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^D ■ The neutral Higgs pair production at the LHC is studied in the MSSM, the 

O 

^^ . large extra dimensional (ADD) model and the Randall-Sundrum (RS) model, 

'j^. where the total cross section can be significantly enhanced compared to that in 

Oh! 

D . the SM. The Pj,, invariant mass and rapidity distributions of each model have 

^ ■ 

been shown to be distinctive: The ADD model raises the px and invariant 

mass distributions at high scales of pt and invariant mass; in the RS model 

resonant peaks appear after the SM contribution dies away; the SM and the 

MSSM distributions drop rapidly at those high scales; in the ADD and the 

RS models the rapidity distributions congregate more around the center. It is 

concluded that various distributions of the Higgs pair production at the LHC 

with restrictive kinematic cuts would provide one of the most robust signals 

for the extra dimensional effects. 

PACS numbers: 04.50.-hh,12.60.-i,13.85.-t,13.90.-hi 



Typeset using REVTeX 



I. INTRODUCTION 

The standard model (SM) has been very successful in explaining experimental signals, 
including recent results from the CERN LEP-II collider |l[]. Nevertheless one of the most 
important ingredients of the SM, the Higgs mechanism, has not been experimentally probed 
yet; it is responsible for spontaneous electroweak symmetry breaking accompanying the mass 
generation of the W^ and Z^ gauge bosons as well as the SM fermions. Recently, ALEPH 
group of the LEP-II has reported the observation of an excess of 3a in the search of the 
SM Higgs boson, which corresponds to the Higgs mass about 114 GeV [0. As operations of 
the LEP-II have been completed, the decision, whether the observations are only the results 
of statistical fluctuations or the first signal of the Higgs boson production, is suspended 
until the Fermilab Tevatron II and/or the CERN LHC experiments running [0]. Thus it is 
naturally anticipated that primary efforts of future collider experiments are to be directed 
toward the search for Higgs bosons 0]. 

In particular at hadron colliders, the pair production of Higgs bosons holds a distinctive 
position in understanding the Higgs mechanism |^. First, it may provide the experimental 
reconstruction of the Higgs potential, as the triple self-coupling of Higgs particles is involved. 
The establishment of the Higgs' role in the electroweak symmetry breaking is crucially 
dependent on this measurement. Second, the signal-to-background ratio is significantly 
improved compared to that of a single Higgs boson production. The invariant mass scale 
of the single Higgs production is fixed by the Higgs mass, of order ~ 100 GeV. Thus their 
detection through heavy quark decay modes suffers from large QCD backgrounds. Besides, 
one viable decay mode /i — > 77 has a very small branching ratio of order 10^'^ 0. For 
the pair production of the Higgs particles, the four 6-jets in the final states are energetic, 
reducing the main background hbb with soft 6-jets [0. Third, this is a rare process in the 
sense that the effects of physics beyond the SM can remarkably enhance the cross section 
with respect to that in the SM; the minimal supersymmetric standard model (MSSM) ^ 
provides some parameter space for the large enhancement of the total cross section, which 



should accommodate the large Yukawa coupling of the b quarks, the resonant decay of 
H —>■ hh, and/or dominantly large contribution of the squark loops [0; extra dimensional 
models provide tree level diagrams mediated by the Kaluza-Klein (KK) gravitons, leading 
to large total cross sections. In fact, these new theoretical approaches have drawn extensive 
attention as candidates for the solution of the gauge hierarchy problem, the existence and 
stability of the enormous hierarchy between the electroweak and Planck scales. Therefore, 
it is worth studying the production of a neutral Higgs pair with the effects of the MSSM and 
the extra dimensional models, and finding the characteristic distribution of each model. We 
shall restrict ourselves to the procedure at the LHC, which is scheduled to start operating in 
2005. Despite the assurance of the detectibihty of the lightest CP-even MSSM Higgs boson 
with y/s > 250 GeV and / Cdt > 10 fb^^ [^], the LHC takes a practical advantage over the 
future e~^e~ linear colliders which lack specific plans yet. 

The remainder of the paper is organized as follows. The next section details the neutral 
Higgs pair production at the LHC in the SM, the MSSM, the large extra dimensional (ADD) 
model [|TU]|, and the Randall-Sundrum (RS) model [[T^]. In Sec. HI, we discuss various 



distributions useful to discriminate effects of each model from others. The last section deals 
with a brief summary and conclusions. 

II. THEORETICAL DISCUSSIONS AND FORMULAE 

The production of a Higgs boson pair at a hadron collider proceeds through several 
modes; WW fusion, bremsstrahlung of Higgs bosons off heavy quarks, and gluon-gluon 
fusion. At the LHC, the gg fusion is expected to play a main role, since the gluon luminosity 
increases with beam energy. In this paper, we focus on the process gg — > hh, where the h is 
the lightest Higgs boson in each theory. 

The invariant amplitude squared is generally written as, in terms of the helicity amplitude 
A^AiA2 for the initial gluon helicities Ai(2) 0, 



i^r = 2 ■ - ■ - • - [|A^++i^ + |A^+_r + i-M-+r + \M-\'\ , (1) 



where the factor 2 refers the color factor [Tr(T'^T'')]^ = 2, the factor 1/4 the initial gluon he- 
licity average, the factor 1/64 the gluon color average, and the final factor 1/2 the symmetry 
factor for the two identical Higgs boson production. For the processes with CP conservation, 

the helicity amplitudes are related as Ai++ = ^A and 7Vl+_ = A^_+. 

For the numerical analysis, we use the leading order MRST parton distribution functions 
(PDF) for a gluon in a proton ||12|. The QCD factorization and renormalization scales Q are 



set to be the hh invariant mass, i.e., Vs. The Q^-dependence is expected to be small on the 
distribution shapes, which are of our main interest. The center-of-momentum (cm.) energy 
at pp collisions is ^/s = 14 TeV. And we have employed the kinematic cuts of Pj, > 25 GeV 
and \ri\ < 2.5 throughout the paper. 

A. In the SM 

In the SM, there are two types of Feynman diagrams as depicted in Fig. 1. One is the 
triangle diagram where a virtual Higgs boson, produced from gg fusion through heavy quark 
triangles, decays into a pair of Higgs bosons. The other is the box diagram where the Higgs 
pair is produced through heavy quark boxes. It is to be noted that the triangle diagram 
incorporates the triple self-coupling of Higgs bosons. For analytic expressions of all the one- 
loop helicity amplitudes of the process gg — > hh, we refer the reader to Refs. [^|| . Fig. 2 
shows the total cross section of the SM Higgs pair production at the LHC, as a function of 
the Higgs mass m/j. At nih — 100 GeV, the atot is of order 60 fb, which decreases rapidly 
with increasing rnh- 

B. In the MSSM 

The existence of a fundamental scalar particle in the SM causes the well-known gauge 
hierarchy problem. Traditional approaches are to introduce new symmetries, motivated by 
the chiral symmetry for light fermion masses and the gauge symmetries for gauge boson 
masses. Supersymmetry is one of the most popular candidates for this new symmetry. 



The MSSM Higgs pair production undergoes distinctive contributions coming from the 
new Higgs and squark sectors, which provide the possibihties to greatly enhance the total 
cross section of the process. First, two doublets and thus two different vacuum expectation 
values of the Higgs fields allow the Yukawa coupling of the b quark to be compatible with 
that of the top quark, which corresponds to the large tan/5 case. Much smaller mass of 
the b quark may makes its loop contribution even larger than the top quark contribution 
p. Second, two Higgs doublets imply the presence of a heavy CP-even neutral Higgs boson 
H, which can decay into two light Higgs bosons if kinematically allowed. This resonant 
contribution gg ^ H —>■ hh is shown to enhance the total cross section with respect to 
the SM case by about an order of magnitude. Third, the MSSM permits a parameter space 
where the 6 or t loop contributions can exceed the SM quark loop contributions by more than 
two orders of magnitude . For the maximization of squark loop contributions which occurs 
through the b loops, this parameter space should allow large value of tan/3, considerably 
light 61 mass, and large mass of A and/or |/i|. Since the third enhancement possibility has 
rather restrictive parameter space (for example, the squark loop contributions are practically 
negligible unless mi < 120 GeV), we consider only the quark loop contributions, as in 
Ref. [^. The corresponding Feynman diagrams are the same as the SM case, except for 
the different coupling strengths and the presence of the neutral heavy Higgs boson H (see 
Fig. 1). 

The MSSM total cross section as a function of the Higgs mass shows similar behavior to 
that of the SM, except for the overall enhancement, as shown in Ref. 0. As anticipated, 
the total cross sections in the large and small tan/? cases are much increased compared to 
the SM case. In particular, the large tan/? value with mh — 100 GeV leads to an order of 
magnitude enhancement of the cross section. 



C. In the ADD model 

The gauge hierarchy problem has been approached without resort to any new symmetry 
by Arkani-Hamed, Dimopoulos, and Dvah (ADD) [|T0|. One prerequisite of the hierarchy 
problem itself is removed; the Planck mass is not fundamental; the nature allows only one 
fundamental mass scale Ms which is at the electroweak scale. By introducing the N > 2 extra 
dimesional compact space, the observed huge Planck mass is attributed to the large volume 
of the extra space, since Mpj ~ Mg "^^i?^ where the R is the size of the extra dimension. The 
prohibition of the SM particles' escaping into the extra space can be achieved by describing 
the matter fields as open strings of which the end-points are fixed to our four- dimensional 
world. 

Of great interest and significance is that this idea is testable at colliders. The KK 
reduction from the whole (4 + A^)-dimensions to our four- dimensional world yields towers of 
massive KK-states in the four-dimensional effective theory. Even though the coupling of a 
KK-graviton to the ordinary matter fields is extremely suppressed by the Planck scale, tiny 
mass-splitting ArrixK ~ 1/-R (which is about 10^^ eV for the N = 2 case) induces summation 
over all KK-states, which compensates for the Planck scale suppression. In addition to single 



graviton emission processes as missing energy events |T3|, the indirect effects of the massive 



graviton exchange on various collider experiments |T^, and the possible Lorentz and CPT 
invariance violations through the change of the metric on the brane |jl5| have been extensively 
studied. 

For the Higgs pair production through the gluon-gluon fusion, there exists a tree level 
Feynman diagram mediated by spin-2 KK-gravitons (see Fig. 3). Based on the effective 



four-dimensional Lagrangian |ITqJT^, the helicity amplitudes are obtained as 



M++ = M— = , (2) 



M+- = M-+ = ^{mt-ut) 



16A, 4 

Mi 



which are to be added to those in the SM. The interference effects between the ADD and 



the SM are proportional to 1/M|, which are, at the energy scale below Ms-, dominant over 
pure ADD effects. 

There might be a concern with our ignorance of the parton mode qq —>■ hh mediated 
by the KK-gravitons. The concern appears reasonable: The characteristic parton energy 
scale Vs of the process in extra dimensional models is of the order Ms ~ TeV, unlike a few 
hundred GeV scale in the SM and the MSSM; the dominant momentum fraction x may not 
be so small as in the SM and the MSSM cases and the magnitude of the parton distribution 
functions of, in particular, the valence quarks becomes substantial. In the following, we 
have taken into account of the parton mode qq —>■ hh, which has the scattering amplitude 
squared as 

1 



iMl^^qq ^ hh) = ^t - u)\tu - ml) . (3) 

According to our numerical analysis, the contribution of qq -^ hh mode turns out to be 
at most a few percent to the total cross sections, because the amphtude squared itself is 
smaller than that of gg -^ hh mode (by a factor of about 1/36) while, at TeV Vs, the PDF 
of a valence quark and a sea quark is of the same order with that of two gluons. 

Note that in the parton cm. frame, the nonzero helicity amplitude can be written as 

•^+-lgg-c.m. = -]^PrS, (4) 

where the p^ is the transverse momentum of an outgoing Higgs particle. Unitarity is appar- 
ently violated at high energies, which is expected from the use of an effective Lagrangian. It 
is reasonable that we only consider the region where our perturbative calculations are reli- 



able, which can be achieved by excluding the region with high invariant mass. In Ref. [jTH 
the partial wave amplitudes of the elastic process 77 — > 77 in the ADD model have been 
examined, yielding the bound on the ratio M^/vs, and a valid region is conservatively 
found to be Vs < O.9M5. In the following analyses, we impose an additional kinematic 
bound such as 

M^^h < 0.9 Ms . 



In Fig. 2, we present the total cross section as a function of rrih with Ms = 2.5 TeV. As 
expected from the presence of a tree level diagram in the ADD model, the total cross section 
is substantially increased with respect to the SM case. For a pair production of Higgs bosons 
with mass 115 GeV, we have obtained 

asM+ADD(mfe = 115 GeV, Ms = 2.5 TeV) _ ^ ^^ (q) 

asuimh = 115 GeV) ^ ' ^ ^ 

We note that the ADD effects lead to a gentle drop-off of the o"tot with respect to the Higgs 
boson mass, contrary to the rapid decrease in the SM. Thus, if the Higgs mass is large, the 
ADD model could still produce substantially large number of Higgs pairs at the LHC unlike 
in the SM case. 

D. In the RS model 

More recently, Randall and Sundrum (RS) have proposed another extra dimensional 
scenario where, without the large volume of the extra dimensions, the hierarchy problem is 



solved by a geometrical exponential factor, called a warp factor |11|. The spacetime in this 



model has a single S^ / Z2 orbifold extra dimension with the metric 

ds^ = e-^''''^^^^r]^^dx^dx'' + rld(l)^, (6) 

where the is confined to < |0| < vr. The Vc is the compactification radius which is to 



be stabilized by an appropriate mechanism [19|. Two orbifold fixed points accommodate 



two three-branes, the hidden brane at = and our visible brane at |0| = vr or vice versa. 
The allocation of our brane at |0| = vr renders a fundamental scale rriQ to appear as the 
four-dimensional physical mass m = e~^^'''^mo, which answers the hierarchy problem. And 
the effective Planck mass is 



M|i = {M^/k){l-e 



~2krciT\ 



where the M is the five-dimensional Planck scale. Note that all of the Mpi, fc, and M are of 
the Planck scale. 
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The compactification of the fifth dimension leads to the following interaction Lagrangian 
in the four- dimensional effective theory |pO||, 



1 1 °° 



where A^ = e~^'^'''" Mp\. Unlike almost continuous KK-graviton spectrum in the ADD model, 
we have one zero mode of the KK-gravitons with the coupling suppressed by the Planck scale, 
and the massive KK-graviton modes with the electroweak scale coupling A^. The masses of 
the KK-gravitons are also at electroweak scale, given by ||21|| , 



k 

Mpi 

where the x„'s are the n-th roots of the Bessel function of order one. 

The scattering amplitudes of the KK-mediated diagrams in the narrow width approxi- 
mation can be derived from the ADD ones with the following replacement in Eq. (2) pO|l : 



A__J_y- I (9) 

where the total decay width of the n-th KK-graviton is r„ = pm„x,^(A;/Mpi)^, and the p, 
fixed to be one, is a model-dependent parameter |2^ . 



The observables based on the four- dimensional effective theory are determined by two 
parameters, (A^r, A;/Mpi). The value of k/Mpi may be theoretically constrained to be less 
than about 0.1 |2^: The magnitude of the five-dimensional curvature, -R5 = —20k^, is 
required to be smaller than M^ (~ M|j), so that the classical RS solution derived from the 
leading order term in the curvature remains reliable. The A^^ is expected to be below 10 TeV 
in order to explain the hierarchy problem. Unlike the Ms in the ADD case, the A^^ does not 
play the role of a cut-off, relieving the concern of the unitarity violation. According to the 
phenomenological studies of the cross section of e^e^ —>■ fi^p^ in the RS model, only the 
case with large value of k/Mpi hints the unitarity violation; even for k/Mpi ~ 1, unitarity 
violation can occur at cm. energy of several TeV [^; the current LEP-II experiments and 
the Tevatron run-I have provided a lower bound of A,r to be about 1.5 TeV in the case of 
k/Mpi = 0.1. 
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In Fig. 2, we plot the total cross sections with respect to the Higgs mass within the RS 
model. We set A^^ = 3 TeV and k/Mp\ = 0.1. For the equity in comparing with the ADD 
case, we have employed the upper bound of Mhh < 0.9 A^. Though smaller than in the ADD 
case, the total cross section in the RS case is larger than that in the SM: 

c^SM+Rs(m;, = 115 GeV, Ms = 3 TeV) 

asM(m, = 115GeV) ^ ' ' ^ ^ 

The rate of the drop-off of atot against rrih is similar to the SM case. 

In order to demonstrate the dependence of k/Mpi and Ayr, Fig. 4 shows the total cross 
section within the RS model as a function of Ayr, considering three values of the ratio 
k/Mp\ = 0.01, 0.1 and 0.3. The Higgs mass is set to be 100 GeV and the upper bound 
in Mhh is not applied. As the A,r increases, the o"tot drops rapidly. And it can be seen 
that smaller value of the ratio k/Mp\ produces larger cross section. This is due to that the 
amplitude squared in the narrow width approximation is inversely proportional to {k/MpiY 
at each resonance, which yields dominant contribution. 

III. NUMERICAL DISCUSSIONS AND DISTRIBUTIONS 

In the previous section, we have shown the possibilities that the Higgs pair production 
can be greatly enhanced at the LHC. In such a circumstance, it is worthwhile to search 
for appropriate distributions which enable us to distinguish the contributions of one model 
from others. In the numerical analysis of the distributions, we have employed the following 
parameters: The Higgs mass is set to be 100 GeV; in the MSSM, tan/5 = 30, fi = — 640 
GeV, Mi = Mi = 1000 GeV, and At = A^ = - 410 GeV; the Ms in the ADD model is 2.5 
TeV; in the RS model, A^ = 3 TeV and k/Mpi = 0.1. 

In Fig. 5, we present the pj-distributions in the SM and the ADD cases. While the 
SM Pj,-distribution peaks around 150 GeV and drops rapidly with increasing Pt, the ADD 
effects slowly draw the distribution up at high pt region. Note that the absence of the 
differential cross section aX px > I TeV is due the employment of the upper bound in the 
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Mhh- Fig. 6 shows the pT-distributions in the SM and the RS cases. The presence of the 
Kaluza-Klein gravitons in the RS model leaves a clear shape of resonance. We caution 
readers that the pr-axis is plotted by log-scale: Even if the da^^~^^^ / dpr at px > 300 GeV, 
which practically vanishes in the SM, is smaller than that at pt < 300 GeV by an order 
of magnitude, the extensive contribution at high pt region renders the total cross section 
substantially enhanced. In Fig. 7, we demonstrate the MSSM pr-distribution for the large 
tan/3 case. It drops rapidly with increasing pt as in the SM case, while it peaks around 
25 GeV, lower than the SM case. The magnitude of the differential cross section is shown 
about twenty times larger than that in the SM case. 

Fig. 8 illustrates the invariant mass distributions of the Higgs pair in the SM and the 
ADD cases. The SM case, where the top quark loop contributions are dominant, peaks 
around the threshold \/S ~ 2?Tij. The ADD effects gently raise the M/i/j-distribution at high 
Mhh region. It is expected that the blind application of the effective Lagrangian in Eq. (H) 
without any upper bound in Mhh would yield continual increase in the M/^/i-distributions, 
which would show an apparent violation of unitarity. We display the M/i/^-distribution of 
the RS case in Fig. 9. The high resonance peak in addition to the SM distribution implies 
the first KK-state of gravitons with rrii ~ 750 GeV. While at the e~^e~ -^ fi^fx^ process 
the KK-gravitons appear as almost regularly spaced peaks |^, the hadronic convolution 



of the parton level processes obscures successive and separated peaks as the classical KK- 
signature. Fig. 10 displays the M^/i-distribution of the MSSM case with large tan/?, which 
undergoes dominant contributions from b quarks. Thus a peak appears just above the 
kinematic threshold, \/§ ~ 2mh. 

Finally, we illustrate the rapidity distributions in Fig. 11. It can be seen that the extra 
dimensional models produce a Higgs pair somewhat more centrally in rapidity than the SM 
and the MSSM do. More restrictive cut on rj, such as r/ < 1.0, would eliminate a substantial 
portion of the SM and the MSSM contributions. 



11 



IV. CONCLUSIONS 

The pair production of neutral Higgs bosons from the gluon-gluon fusion at the LHC 
has been studied in the SM, the MSSM, the large extra dimensional model and the Randall- 
Sundrum model. We have shown that both the supersymmetric and extra-dimensional 
models can substantially enhance the total cross section of the Higgs pair production. In 
the MSSM case, the large tan/3 value makes possible the 6-quark contribution dominant over 
the top quark contribution and the resonant decay of a heavy Higgs particle into two light 
Higgs particle. The extra dimensional models allow tree level Feynman diagrams mediated 
by the Kaluza-Klein gravitons, which significantly increase the total cross section. Since 
the ADD model has been shown to undergo the violation of partial wave unitarity at high 
energies, we have employed an upper bound in the invariant mass of two Higgs bosons, 
which is obtained from the analysis of the J-partial wave amplitudes of the elastic process 
rjrj _^ rjrj_ \t^ addltlou. It was shown that the total cross section in the ADD case with the 
upper bound in Mhh decreases gently with increasing Higgs mass, whereas those in the SM, 
the MSSM and the RS cases decrease rapidly. 

If Higgs pairs are produced at hadron colliders much more than the SM prediction, the 
three non-standard models considered here are good candidates for the explanation. We 
have demonstrated the p^, invariant mass and rapidity distributions of each case. The 
distribution shapes are shown to be different for each model, providing valuable criteria to 
distinguish the contribution of one model from others. The p^-distribution in the SM peaks 
and drops rapidly; the MSSM one for the large tan/9 case peaks just above the threshold 
of pt and also drops rapidly; the ADD effects induce a slow raising after the SM peak; 
the RS effects can be discovered by the presence of a resonance peak. The invariant mass 
distributions are similar to the j9j,-distributions: The SM and the MSSM have peaks around 
a few hundred GeV; the ADD effects gently raise the M/j/^-distribution at high energies; the 
RS contribution yields a series of resonant peaks. Therefore, restrictive cuts on the pt and 
Mhh would eliminate the main contributions of the SM and the MSSM cases, which provides 
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one of the most straightforward methods to signal the existence of low scale quantum gravity 
effects. Finally, the rapidity distributions in the ADD and the RS models show significantly 
narrow peaks around 1] = 0, which implies the large contributions at high p^ region. 
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FIG. 1. The Feynman diagrams of the gg -^ hh process in the SM (MSSM). 
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FIG. 2. The total cross section of the Higgs pair production as a function of the Higgs mass 
at the LHC with -^s = 14 TeV in the SM, the ADD and the RS cases. The string scales are set to 
be Ms = 2.5 TeV and Ayr = 3 TeV for the ADD and the RS cases, respectively. The upper bounds 
in Mhh < 0.9Ms{At,) are employed. 
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FIG. 3. The Feynman diagrams of the gg -^ hh in the ADD and the RS models. 
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FIG. 4. The total cross sections only with the RS effects as a function of Ayr, considering three 
values of the ratio k/Mpi = 0.01, 0.1 and 0.3. 
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FIG. 5. The Pj. distributions of the Higgs pair production in the SM, and the ADD cases. 
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FIG. 6. The Pj, distributions of the Higgs pair production in the SM and the RS cases. 
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FIG. 7. The Pj, distributions of the Higgs pair production in the SM, and the MSSM. 
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FIG. 8. The M/j/j distributions of the Higgs pair production in the SM and the ADD cases. 
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FIG. 9. The M^h distributions of the Higgs pair production in the SM and the RS cases. 
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FIG. 10. The M^^ distributions of the Higgs pair production in the SM and the MSSM. 
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FIG. 11. The ri distributions in the SM, the MSSM, the RS and the ADD cases. 
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